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Description 

Design Technique for Use in 
Engineering of Deamidation Rates of 
Peptides, Proteins, Hormones, and 
Peptide-Like, Protein-Like and 
Hormone-Like Molecules 

Background of Invention 

[0001] jhe deamidation of peptides and proteins as well as 

molecules related to peptides and proteins is a well known 
phenomenon. In this reaction Asn or Gin residues are 
gradually changed into Asp and Clu residues respectively 
as well as their isomers and charged products. The rate of 
this reaction is dependent on the primary sequence, 
three-dimensional structure, pH, temperature, buffer 
type, ionic strength and other solution properties. The 
half-time varies from less than 1 day to more than a cen- 
tury. The reaction introduces a negative charge into the 
molecule. In addition, the isomerization products fc-Asp 



and fc-Glu as well D-isomerized forms and chain cleavage 

also accompany the reaction. 

[0002] The stability of Asn and Gin in pharmaceutical and other 

types of commercial preparations is a major field of study. 

Efforts have been made to discover formulation conditions 

that will control the rate of deamidation of amides in 

these preparations. 
Summary of Invention 

[0003] The inventions described here allow the engineering of 
molecules with specific amide structures that will deami- 
date at a specified rate. This procedure can be used to 
design stable and unstable forms for pharmaceutical, in- 
dustrial, and other products. This can be used to increase 
the shelf-life of such products through minor modifica- 
tions, prevent or at least slow down the gradual formation 
of impurities in preparations with these modifications, and 
may make possible as a result of minor modifications the 
use of products that would otherwise be too unstable for 
practical purposes. The engineering of products with un- 
stable amides that are programmed to deamidate at spe- 
cific rates may also be a valuable application of this pro- 
cedure. 

[0004] This invention is based on newly discovered rules that 



govern the deamidation rates of amides. These rules can 
be derived using a newly discovered method. Many of 
them have been derived and more can be derived using 
this method. These rules involve the effects of steric and 
catalytic factors in the nearby proximity of amides. It has 
been shown that these rules can be quantitatively applied 
to determine the deamidation rates of amides to a very 
high accuracy. Even without the exact use of these rules, 
the qualitative effect on the deamidation rate of adding or 
subtracting various functional groups in various positions 
has been discovered. 

[0005] This procedure makes available a wide variety of possible 
variations that will modify deamidation rates. Modification 
of the residue to the right of the amide using either natu- 
ral or non-natural amino acids, as well as modification of 
Asn or Gin and determination of the quantitative effects 
on the deamidation rate are now available. 

[0006] insulin is an example of an important pharmaceutical 

molecule that deamidates during storage. The procedure 
described in this patent allows the invention of stable 
forms of insulin that require a minimum of modifications 
and will likely have the same biological activity as the na- 
tive form. Two specific modifications are given as exam- 



pies. Many other pharmaceutical preparations, such as in- 
terferon and tumor necrosis factor are thought to be de- 
activated by deamidation. 
Detailed Description 

[0007] General Method: 

[0008] it W as known before the invention of the method de- 
scribed here, that the sequence and structure around an 
amide has a large effect on the deamidation rate. What 
was unknown was how this effect worked or any quantita- 
tive information that would allow engineering of stable or 
unstable amides or amides with predetermined rates. As a 
result of experiments that I carried out, sufficient data to 
allow such a method to be invented were collected. The 
method described here uses this data and makes clear 
how more such data can be applied to further extend the 
method. It has been developed using Asn and Gin 
residues in peptides. It applies to peptides, hormones, 
proteins, and peptide-like, hormone-like, protein-like 
molecules, as well as other molecules that have unstable 
amides in similar surroundings. It also applies to isomer- 
ization and racemization of Asp and Glu in similar circum- 
stances to Asn and Gin. No such system has, to my 



knowledge, ever been described. 

[0009] This method involves two parts. The first of these is the 
use of experimental data to determine the effect of atoms 
and structures made up of atoms on the effective activa- 
tion energy of the deamidation reaction. Each type of 
structure is assumed to increase or decrease the overall 
activation energy of the reaction and is assigned a coeffi- 
cient depending on its location in the structure. The sec- 
ond part involves the use of these derived coefficients to 
predict primary sequence deamidation rates for use in en- 
gineering amides with particular rates or stable structures. 
Even without the quantitative basis, this method makes 
clear the required modifications necessary to achieve a 
particular qualitative change in a deamidation, isomeriza- 
tion, or racemization rate. 

[0010] t wo types of deamidation are obviously present. The one 
on which this method is based, and which is most preva- 
lent for amides with half-times less than a few hundred 
days, depending on conditions and providing especially 
catalytic ions are not present, is most strongly effected by 
the structure to the right of the amide (e.g. in the se- 
quence GlyXxx(Amide/Acid)YyyCly the identity of Yyy is 
the most important factor). Also present is at least one 



more mechanism that is usually slower and has different 
sequence dependence. It is possible that this dependence 
as well as the left hand structure dependence (Xxx in the 
sequence GlyXxx(Amide/Acid)YyyGly) can also be modeled 
with a similar system, but this has not yet been demon- 
strated. 

[001 1] M 0S t of the structures that increase the activation energy 
and therefore slow the rate, apparently do so through 
steric hindrance of the reaction. Sometimes both a cat- 
alytic and a steric effect are present in a particular atom or 
group of atoms. 

[0012] This model is based on the understanding that an individ- 
ual molecule deamidates only as a statistical phe- 
nomenon. It does not deamidate at a particular rate. Each 
amide has a probability that it will deamidate over a par- 
ticular time interval. The combined statistics of an ensem- 
ble of such molecules produce the deamidation rate. For 
example, the deamidation rate of ClyXxxAsnGlyGly is 
about 1 day, so each individual molecule has a 50% 
chance of deamidating in 1 day or in one second only a 
0.00058% chance of deamidating. The deamidation of 
even a fast amide occurs infrequently. 

[0013] Regardless of the details of the mechanism, this method 



assumes the amide must adopt a particular conformation 
in order to react. Sterically hindering groups have at least 
two ways of interfering. Any atom or group of atoms in- 
terfering at the reaction site can increase the energy re- 
quired for the amide to move into position and raise the 
activation energy. Also, when such a structure can adopt 
more than one conformation it may interfere to a signifi- 
cant extent only a certain fraction of the time. Steric ef- 
fects, therefore, depend on group size, the conformations 
available to the group, and the distribution function of oc- 
cupation of these conformations. These ideas are used to 
understand steric effects. 
[0014] Calculation of Constants: 

[0015] Median values for carboxyl side residues of Asn and Gin 

were used in these derivations. The measurements for Asn 
are the most reliable and are used for these derivations. 
Gin calculations are likely to contain significantly more 
experimental error, as well as the possibility of contribu- 
tions from other types of deamidation mechanisms, but 
are still commercially useable. An average correction for 
the effect thought to be the largest, hydrolysis, has been 
made to all rates, but only significantly affects the Gin 
rates and very slow Asn rates. The half-time of hydrolysis 



was assumed to be 8010 days, but may be faster or 
slower. Regardless of the derivation based on Asn, calcu- 
lations of Gin based on these values are found to approxi- 
mately agree with the currently available data. 
[0016] The natural log of the median hydrolysis-corrected rate 
constants times 100 (s -1 ) for the Asn and Gin pentapep- 
tides as a function of carboxyl-side residue as computed 
from the experimental data are listed in Table 1. The first 
columns are (100)(ln(k)) and are directly proportional to 
the Arrhenius activation energy. The second columns are 
normalized to Gly. This leaves only the side-chain effects 
and also defines the effect of the Gly hydrogen side-chain 
or a similar atom with mostly steric effects in this p- 
position as zero. 



Table 1. (100)ln(k) deamidation for the medians of Xxx in GlyXxx(Asn/Gln)YyyGl> 



Asn s 

Calculated Values Normalized to Gl\ 



Calculated Values Normalized to Gly 



Gly 
His 
Ser 
Ala 
Asp 
AmCys 
Thr 
Cys 
Lys 
Met 
Glu 
Arg 
Phe 
Tyr 
Trp 
Leu 
Val 
He 



-1186.3 
-1405.8 
-1448.3 
-1493.1 
-1520.1 
-1528.4 
-1556.9 
-1566.2 
-1579.4 
-1582.8 
-1587.4 
-1587.0 
-1598.3 
-1611.5 
-1630.3 
-1652.6 
-1724.8 
-1747.4 



0 

219.4 

262.0 
306.8 
333.7 
342.1 
370.6 
379.8 
393.1 
396.5 
401.1 
400.7 
411.9 
425.1 
444.0 
466.3 
538.5 
561.1 



-1831.1 
-2181.5 
-2165.6 
-2178.2 
-2393.4 

-2136.5 
-1958.8 
-2184.9 
-2104.7 
-2313.2 
-2290.7 
-2433.6 



-2198.9 
-2230.9 
-2210.2 



0 

350.4 
334.4 
347.1 
562.3 

305.3 
127.6 
353.7 
273.5 
482.0 
459.5 
602.5 



367.7 
399.7 
379.0 



§ A correction of 8,0 1 0 = half-time for hydrolysis was applied to the experimental values before 
normalization to Gly. # k = sec" 1 . 



[001 7] E ff ect of Aliphatic Side-Chains: 

[0018] Aliphatic side-chains are assumed to have only steric ef- 
fects. Any catalytic effects that might be present are 
masked by the steric effects. The calculations of constants 
for groups in each position are done by finding the differ- 
ence in activation energy on structures for which a group 
was substituted. For the purposes of this report the a- 
carbon is the backbone carbon in the peptide chain with 
successive positions assigned Creek letters as shown in 
Figure 1. Figure 2 illustrates these calculations for three 
different substitutions. In some cases more than one 
method of calculation is possible. In all such cases the 



values are found to be similar and are averaged. Continu- 
ing with the method shown in Figure 2 we can make the 
following calculations: 
[0019] (Ala Gly) = 306.8 0 = 306.8 

[0020] Y -CH 3 : (Val Ala)/2 = (538.5-306.8)/2 = 115.9 

[0021] Y -CH 3 : (Thr Ser) = 370.6 - 262.0 = 108.6 

[0022] Average y-CH^ (115.9 + 108.6)/2 = 112.3 

[0023] 8_CH 3 : (lie Val) = 561.1-538.5 = 22.6 

[0024] 6_CH 3 : (Leu (Ala + y-CH^) = 466.3 306.8 112.3 = 23.6 

[0025] Average 5-ChM (22.6 + 23.6)/2 = 23.1 

[0026] Substitution increments for replacement of a hydrogen 
atom by a methyl group in the 3, y, and 5 positions are 
therefore 306.8, 112.3, and 23.1. The internal consistency 
with two redundant examples verifies the procedure. The 
values become successively smaller as the methyl groups 
are placed farther from the reaction site as expected. 

[0027] These values can also be derived from the statistical 

probability that the methyl groups will be in a position to 
interfere with deamidation. Define the effect of the B-CH 

K 3 

group as X, and substitute one of its hydrogen atoms with 
a y-CH group. Reference to a CPK model shows that the 



new methyl group will occupy a similarly interfering posi- 
tion for about 120° out of the 360° rotation around the a- 
3 bond. Thus the effect of the second methyl group 
should be about (1/3)(X). Examination of CPK models with 
the addition of a &-CH 3 group indicates that the interfer- 
ence will be the same as that of the Y-CH 3 except that it 
occurs for 180° out of the 360° of the tx-p bond rotation, 
but only during 120° of the 360°p-y rotation. Thus the ef- 
fect of replacing a hydrogen atom with a &-CH 3 group is 
(180/360 - 120/360)(120/360)X = (3/6 - 2/6)(l/3) = 
(1/18)X. 

[0028] since X is 306.8, the effects of the y and 8 methyl groups 
should be 102.3, and 17.0 respectively. In fact they are 
112.3 and 23.1. Now the hydrogen atom and carbon atom 
effects can be separated. For example, the substitution of 
a y-CH 3 group involves the removal of a y-H and replace- 
ment by a y-C and 3 5-H atoms. The individual effects 
can be separated. The experimental results can be encap- 
sulated in four equations: 

[0029] h' 3 = 0 (By definition from normalizing to Gly) 

[0030] C + 3 H - H = 306.8 

3 Y 3 

[0031] C + 3 H - H = 112.3 

Y 5 y 



[0032] C + 3 H - H = 23.1 

6 g 8 

[0033] jo a jd j n solving these equations the carbon atom effect is 
postulated to be roughly a constant, "a", times the effect 
of a hydrogen atom if they are in identical positions. While 
not a perfect approximation it is sufficient for most pur- 
poses. Thus: 

[0034] c = a(H ) 
3 3 

[0035] c = a(H ) 

Y Y 

[0036] c = a(H ) 
8 8 

[0037] c = a(H ) 

e e 

[0038] jo determine the value of "a", the cube root of the carbon 
atom result is plotted against the substituent number in 
Figure 3b. By repetitive solving of the equations it is found 
that the graph becomes perfectly linear when a = 1.038. 
Using this value for "a": 

[0039] H = 0 
3 

[0040] c = 0 
3 

[0041] h = (306.8)/3 = 102.3 

Y 

[0042] c = a (306.8)/3 = 106.1 

Y 

[0043] h = (H2.3 + (l-a)(306.8)/3)/3 = 36.1 



[OO 44 ] C = a(112.3 + (l-a)(306.8)/3)/3 = 37.5 

5 

[0045] H = (23.1 + (l-a)(112.3 + (l-a)(306.8)/3)/3)/3 = 7.2 

G 

[0046] C = a(23.1 + (l-a)(112.3 + (l-a)(306.8)/3)/3)/3 = 7.5 

G 

[0047] This calculation also extracts the e hydrogen and carbon 
atom effects even though they have not yet been mea- 
sured exactly. Extrapolation of the line from the first three 
points in Figure 3b also generates a value for the fourth 
point at 

[0048] The assumption on which Figure 3b is based is that the 
fractional volume occupied by an atom or group should 
decrease with approximately the cube of the substituent 
number. The unmodified experimental values derived pre- 
viously from only differences in sequence rates are also of 
interest, but are only known in terms of removal of a hy- 
drogen atom and substitution of a methyl group. This plot 
is shown in Figure 3a and has a correlation coefficient of 
0.99996. The first three values are the experimental val- 
ues derived previously, and the last one is computed from 
the results of the calculation. The position of this fourth 
value optimizes at the same value of "a" where the line in 
Figure 3b becomes a perfect fit. 

[0049] Aromatic Rings: The amino acids Phe, Tyr, and Trp contain 



aromatic rings. Figure 4 shows that these rings are con- 
nected in the y position for Phe and Trp. Tyr is similar and 
only differs from Phe by the addition of an OH group. Us- 
ing the previous results the effect of these groups can be 
isolated: 

[0050] Y C H =Phe-(C + 2 H ) = 411.9 (0 + 2(102.3)) = 207.3 

1 5 5 3 Y 

[0051] Y C 5 H 5 OH = Tyr - (C + 2 = 425.1 (0 + 2(102.3)) = 
220.5 

[0052] Y C H N = Trp - (C + 2 H ) = 444.0 (0 + 2(102.3)) = 

8 6 3 Y 

239.4 

[0053] it also turns out that the aliphatic prediction scheme can 
give very accurate results for these aromatic rings without 
measuring them directly. 

[0054] Looking at the composition of the Phe ring in Figure 4, it 
is found that the y C and 5-C are about as free for these 
purposes as in a non-aromatic ring. The e-H is not as free 
to move, but its direction allows it to interfere with 
deamidation. The remaining atoms are directed away from 
the reaction site. Therefore: 

[0055] Y -Phenyl = C +2C + 2H = 106.1 + (2)(37.5) + 

Y 5 e 

(2)(7.2) = 195.5 
[0056] |f the c carbon atoms are included, y-Phe = 210.5 



[0057] Either result agrees closely with the experimental value for 
Y-Phe of 207.3. 

[0058] For Trp the extra effect of the second ring can be added 
to the experimental value of 207.3 for Phe. Specifically it 
is necessary to determine the effect of the T^-H atom. 
When averaged over the many conformations available to 
an aliphatic the value is very small, but this atom in Trp is 
held in one position that is directed toward the reaction 
site. Figure 4 shows that the ^-H atom occupies a position 
relative to the reaction site that is very similar to the 
nearby 8-C. Using the value for a 8-H atom and correcting 
the e-H to a e-C: 

[0059] y-lndole = y-Phenyl + H +C -H = 207.3+ 36.1 + 7.5 

5 e e 

7.2 = 243.7 

[0060] Direct calculation without reference to Phe yields: 

[0061] y-lndole = C + 2C+H+C+H = 106.1 + (2)(37.5) 

Y 5 e e 6 

+ 7.2 + 7.5 + 36.1 = 231.9 

[0062] B 0 th results are in good agreement with the experimental 
value of 239.4. 

[0063] T yr j S identical to Phe except for the OH group in the para 
position on the phenyl ring. The steric hindrance should 
be similar to that of Phe. Tyr turns out to be slightly 
higher at 220.5 vs. 207.3. Why this slight modification in 



rate occurs is unknown, but it is exerting some type of 
catalytic effect, which this calculation has quantitated. 
[0064] Thus the aliphatic side chains are reduced to individual 

atomic contributions. They can be summed to predict the 

3 

deamidation rate for any sp hydrocarbon chain. These 
atomic coefficients also work well for aromatic rings. 
Therefore, the imide deamidation rates for any aliphatic 
Asn amino-side residue can be accurately computed. This 
ability is extremely valuable in terms of its practical quan- 
titative prediction and engineering possibilities and in the 
understanding that will allow qualitative modifications to 
be made and understood. 
[0065] Other Atoms and Charges: 

[0066] other commonly occurring atoms in the ordinary amino 
acid residue side chains include oxygen, nitrogen, and 
sulphur. Many other possible atoms exist for non-natural 
side-chains. In addition to steric hindrance, these atoms 
can have additional characteristics such as hydrogen 
bonding capabilities and positive and negative charges 
that can modify their contributions. 

[0067] with only data from the naturally occurring amino acids 
some of these effects can be described and utilized. Fur- 
ther experimentation and analysis using these procedures 



can lead to understanding the effects of any atom, group 
of atoms, or similar structure both on the immediate side 
chain of the right hand residue and in proximity from 
other effects. It may be found that under certain circum- 
stances the coefficients for particular combinations of 
groups cannot be directly added to determine the rate. 
This will entail a very minor modification of the proce- 
dures. 

[0068] of particular interest is positive charge in Arg, Lys and 

His. Assuming that the steric hindrance of N is similar to C 
we calculate: 

[° 069 ] Catalytic His = His Phe = 219.4 411.9 = -192.5 

[° 07 °] Catalytic Lys = Lys (C + 2H + C + 2H +C + 2H + 

(3 Y Y 8 6 G 

C + N + All others zero) = 393.1 [0 + (2)(102. 3) + 

G £ 

106.1 + (2)(36.1) + 37.5 + (2)(7.2) + 7.5 + 0.2] = - 49.4 
[° 071 ] Catalytic Arg = Arg (C + 2H + C + 2H + N + H + C 

7 3 3 P YY 66gg 

+ 2 + All others zero) = 400.7 [0 + (2)(102.3) + 106.1 
+ (2)(36.1) + 37.5 + (7.2) + 7.5 +(2)(0.2) ] = - 34.6 

[0072] Figure 5 illustrates these catalytic effects as a function of 
the distance from the reaction site. Two other points de- 
rived from rates on peptides with His one and two 
residues removed from the right hand side are included. 

[0073] The y-OH groups of Ser and Thr have coefficients of 55.6 



(Ser is 57.4 and Thr is 53.8, this is the average). Assuming 
a steric effect similar to that of CH, the catalysis of this 
group is 55.6 - 142.2 = -86.6. 
[0074] cv S and Met contain sulfur. In the case of Cys the net ef- 
fect of the SH group is 201.5. Removal of the steric effect 
by assuming it is similar to a CH 2 yields 201.5 178.4 - 
36.1 = -13. For AmCys (hydrogen atom replaced by ac- 
etamidomethyl group) the sulfur atom has a net effect of 

84.4 (assuming no special effects from the ac- 
etamidomethyl group). Correcting for the steric effect 
gives 84.4 - 178.4 = -94. For Met the net effect of S is 
5.5. Correcting for its steric effect yields 5.5 - 52.0 = - 
46.5. 

[0075] The carboxylic acids, Asp and Clu, both have negative 

charges as well as hydrogen bonding capabilities. The net 
effect of the carboxyl is 129.2 for Asp and 15.4 for Glu. 
Correcting for the steric effects results in a catalytic effect 
of 129.2 - 106.1 - 2(37.5) = -51.9 for Asp and 15.4 - 

37.5 - 2(7.5) = -37.1. 

[0076] The derived increments are listed in Table 2. For conve- 
nience in protein engineering, all values listed are the net 
effects of the group without separation of steric and cat- 
alytic components. The catalytic effect can be isolated by 



subtracting the corresponding steric effects. Values listed 
in bold are computed exactly from experimental rates. 
The others are computed based on certain assumptions. 
The bases for the aromatics are the hydrocarbon effects 
shown previously to be quite accurate. The O values are 
based on both carboxylic acids and hydroxyl groups. 
[0077] Qi n Deamidation: 

[0078] This model for Asn deamidation is also applicable to Gin 
deamidation, although Gin deamidation rates have not 
been as well characterized. While there may be important 
differences in the functional group effects for the longer 
Gin chain, the steric effects are similar as illustrated in 
Figure 6. Figure 6a, shows the Gin values without correc- 
tion for hydrolysis vs. the Asn values with hydrolysis re- 
moved, as listed in Table 1. The fundamental difference 
between Asn and Gin rates has been removed by sub- 
tracting the Gin values from GlyXxxGlnGlyGly and the Asn 
values from GlyXxxAsnGlyGly. 

[0079] The difference between GlyXxxGlnGlyGly and GlyXxxAs- 
nGlyGly is 644.8. In a similar manner to the effects of 
steric hindrance the extra -CH 2 - decreases the probability 
of the Gin side chain reacting. One way to look at this is 
by determining the extra entropy associated with the 



longer Gin chain. Based on the standard entropy of for- 
mation of hydrocarbons, the entropy associated with 
adding an extra -CH^- is 40J/(K mol), which corresponds 
to about A100(ln k) = 480 on our scale. This is 74% of the 
actual experimental value. 

[0080] with an average effect of hydrolysis removed, the points 
scatter around the theoretical 1:1 line as shown in Figure 
6b. Without this correction Figure 6a shows the correla- 
tion to be poor. The two dimensional medians of these 
values are shown. 

[008 1 ] u se f n Molecular Engineering: 

[0082] Table 2 is a listing of the currently determined coeffi- 
cients. While the addition of more values to this table 
would be useful, this table already allows prediction of 
deamidation rates based on the structure around the 
amide and protein engineering of many amide rates. 



Table 2. A (100)ln(k) coefficients for calculating deamidation rates*. 




P 


Y 


5 


£ 




H 


-H 


0 


102.3 


36.1 


7.2 


0.18 


0 


-CH 3 


306.8 


214.5 


59.2 


8.1 


0.19 


0 


-CH 2 - 


204.5 


178.4 


52.0 


7.9 


0.19 


0 


-CH — 


102.3 


142.2 


44.7 


7.7 


0.19 


0 


— C — 


0 


106.1 


37.5 


7.5 


0.19 


0 


-C 5 H 5 


284.5 


207.4 


52.9 


7.9 


0 


0 


-C 5 H 5 OH 


— 


220.6 


— 


— 


— 


— 


- QHeN (Indole) 


390.6 


239.4 


60.4 


8.1 


0 


0 


-C 3 H 3 N 2 + (Imidazole) 


— 


14.9 


— 


— 


— 


— 


-S- 





84.4 


5.5 











-SH 


— 


201.5 


— 


— 


— 


— 


-O- 


— 


19.5 


11.5 


-9.6 


— 


— 


-OH 


— 


55.6 


— 


— 


— 


— 


-co 2 


— 


129.2 


18.2 


— 


— 


— 


— mt 








-136.0 





-49.7 


-42.1 


-N 3 CH 5 + (Guanidino) 


— 


— 


— 


-34.2 


— 


— 


-NH 3 + 


— 


— 


— 


— 


-49.7 


— 


*Bold-face values based directly on experimental rates and k = 


sec" 1 . 







[0083] Various groups are listed in addition to individual atoms 
for convenience. To use this table sum up the coefficients 
for the carboxyl-side residue or for non-peptide like 
structures the groups in similar relation to the amide. 
Multiply bonded groups should be treated as described 
for Phe and Trp. The deamidation half-time is computed 
as follows: 

[0084] Asn Peptides t = [(ln(2))/86400]e (( Sum/100) + 1063) 



[0085] cin Peptides t = [(ln(2))/86400]e (( Sum/100) + 18 311) 

[0086] For most Asn peptides the computed result should be 
within 10% or less of the experimental value in 37°C, pH 
7.4, 0.15 M Tris-HCI buffer or similar conditions. Gin 
rates are more affected by hydrolysis and cannot be veri- 
fied to as high an accuracy with the current data. For Gin 
peptides and long lived Asn peptides, correction for hy- 
drolysis should be made. With a hydrolysis value of 8010 
days: 

[0087] t = 1/(1/8010 + 1/t ) 

l/2total l/2calculated 

[0088] As an example take the sequence GlyXxxAsnThrGly. To 
calculate the deamidation half-time from Table 2, it is 
necessary to consider the Thr side-chain. The structure is 
CH(CH 3 )OH. Adding a p-CH, a y-CH^ and a y-OH and 
computing for Asn: 

[0089] AsnThr succinimide t = [(ln(2))/86400]e ((102 - 3 + 214 5 + 55 - 6)/10 ° + 
= 47.2 days 

[0090] including hydrolysis: 

[0091] AsnThr total t = l/((l/47.2) + (1/8010)) = 46.9 days 

[0092] The experimental value is 46.2 days (1). 

[0093] For the non-natural amino acid norleucine in the se- 



quence AsnNle, the side chain is -CH CH CH CH . Adding 

-i 2 2 2 3 3 

the appropriate groups and computing gives: 

mnCMI a mi j *. r/i n\\ mr Ar\r\-i (( 204.5 + 178.4 + 52.0 + 8.1)/: 

[0094J AsnNle succmimide t = [(ln(2))/86400]e 

+ 11.863) nr r , 

= 95.6 days 
[0095] Correcting for hydrolysis: 

[0096] AsnNle total t = l/((l/95.5) + (1/8010)) = 94.4 days 

[0097] Thus the deamidation rate of AsnNle has been determined 
without actually measuring it. 

[0098] For molecular engineering purposes it is often easiest to 
work from an existing structure. Without knowing the fac- 
tors that contribute to the original deamidation rate, the 
effects of a modification can be computed or estimated. 

[0099] As an example, consider the sequence Phe- 

ValAsn(B3)GlnHis of insulin. This peptide deamidates at 
Asn(B3) with a half-time of 136 days at 37°C, pH 7.4 
sodium acetate, 0.1% methyl paraben, 0.7% NaCI. 

[0100] Deamidation of insulin at this position in pharmaceutical 
insulin preparations is a substantial problem. Many at- 
tempts have been made to modify insulin storage condi- 
tions to suppress this deamidation. 

[0101] This computation method provides a procedure for slow- 
ing insulin deamidation with a minimum modification of 



its structure. The Gin side chain is CH CH CONH . Addi- 

22 2 

tion of two methyl groups would change this to CCCH^CH 
CONH . 

2 2 

[0102] without reference to the primary sequence rate of AsnGIn 
or insulin secondary and tertiary structure effects, the 
half-time can be computed from the current experimental 
rate of 136 days. Converting the 136 value to ln(k) and 
subtracting from the AsnGly value yields: 

[0103] -H86.3 - 100[ln[((ln2)/136)(l/60)(l/60)(l/24)] = 478.3 

[0104] The modification involves removing two hydrogen atoms 
and replacing them by two methyl groups. The sum is: 

[0105] 478.3 2 y-H + 2 Y-CH 3 = 478.3 (2)(102.3) + (2)(214.5) = 
702.7 

[0106] Substituting into the deamidation half-time calculation 
yields: 

[0107] t = [(|n(2))/86400]e ( 702 - 7/100 + 11 863) = i 28 3 days 

[0108] This modification should have a half-time increased by 

nearly a factor of 10 with a minimal modification in struc- 
ture and involves no changes in side chain lengths or 
functional groups. This insulin is likely to have biological 
activity equivalent to the rapidly deamidating form. 

[0109] Conclusions: 



[0110] The methods and results described here provide a new 

quantitative understanding of primary sequence control of 
deamidation. With this quantitative understanding the ef- 
fects of modifications in structure around an amide be- 
come clear, and even qualitative modifications can be un- 
dertaken with confidence and without precise calculation. 



